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^  INTERIM  TECHNICAL  REPORT 

INTRODUCTION 

I 

This  is  the  first  interim  report  under  Contract  NoT  F  30602-71-0132 
entitled  "Investigation  of  Laser  Propagation  Phenomena."  This  effort 
is  aimed  at  providing  theoretical  support  to  the  RADC  Laser  Propagation 
Program.  The  report  covers  the  period  January  1*  1971  to  July  1,  1971. 

The  theoretical  support  is  in  the  area  of  linear  atmospheric 
propagation  phenomena  and  microturbulence  statistics.  Other  areas  include 
theoretical  support  to  the  performance  of  propagation  experiments  and 
In  the  interpretation  and  processing  of  the  data  to  ensure  proper  match 
between  theory  and  experiment. 

During  the  past  six  months  effort  has  been  concentrated  in  two 
areas;  examination  of  proper  averaging  times  required  for  propagation 
experiments  and  theoretical  backup  for  phase  structure  function  measure¬ 
ments.  A  bibliography  on  optical  propagation  prepared  previously  was 
also  updated.  The  work  on  averaging  times  Is  in  direct  support  to 
the  experimental  program  at  RADC  because  it  is  intended  to  supply 
information  to  be  used  In  the  data  taking  procedure.  Some  suggestions 
have  already  been  incorporated  into  the  data-taking  procedure.  The 
phase  structure  function  work  consists  of  the  calculation  of  curves  of 
interest  in  experiment  design  and  interpretation.  Specifically  the 
curves  predict  the  contribution  of  various  portions  of  the  path  to  the 
final  measured  value  of  the  phase  structure  function  for  a  spherical 
wave  and  for  horizontal  and  slant  paths  in  the  turbulent  boundary  layer. 
These  will  be  useful  in  Interpretation  of  results  and  In  future  design 
of  experiments.  The  bibliography  Is  an  extension  of  a  previous  work. 

These  items  will  now  be  considered  in  detail  starting  with  the 
work  on  averaging  times! 


AVERAGING  TIMES 

In  any  type  of  measurement  involving  random  quantities,  the 
quantities  of  Interest  must  be  average  values  because  it  is  only  the 
average  values  which  have  any  chance  of  being  reproducible  and  predictable. 
Such  quantities  as  means,  variances  and  covariances,  which  are  of  interest 
in  one  application  or  another  are  all  examples.  However,  in  order  for 
these  quantities  to  be  reproducible,  and  therefore  meaningful,  the 
averages  must  be  taken  over  a  representative  set  of  values,  a  set  large 
enough  to  include  all  the  values  to  be  encountered  in  the  proper  pro¬ 
portion.  It  is  to  answer  the  question  of  how  much  data  to  take  in 
order  to  provide  such  a  set  of  values  that  this  study  has  been  under¬ 
taken.  This  study  is  by  no  means  complete.  The  chief  result  to  date 
has  been  an  examination  of  the  literature  to  determine  the  various 
approaches  used  there.  Results  of  this  examination  will  be  presented 
here. 


1 


There  are  several  criteria  for  determining  what  a  proper 
averaging  time  might  be.  They  are  based  on  precision  desired  and  on 
the  duplication  of  conditions  required  by  assumptions  in  theoretical 
developments.  For  example,  in  theories  of  light  beam  propagation  in 
random  media  the  assumption  of  homogeneity  and  isotropy  are  almost  always 
invoked.  Further  one  may  well  use  the  central  limit  theorem  to  simplify 
predictions  involving  not  too  long  propagation  paths.  In  taking  the 
data  one  generally  uses  the  ergo  die  theorem  and  measures  time  averages 
rather  than  ensenble  averages.  Finally  the  upper  and  lower  spectral 
limits  of  the  desired  data  must  be  considered.  All  of  these  points 
depend  on  the  averaging  time  for  satisfactory  physical  implementation. 

In  this  section  we  present  a  review  of  the  various  approaches  in  the 
literature  to  ques cions  involving  averaging  times,  indicating  assump¬ 
tions  and  making  a  few  limited  observations.  Much  of  the  work  described 
comes  from  the  micrometeorological  literature,  the  rest  stemming  from 
conmuni cation  sources.  The  various  authors  consider  various  aspects 
of  averaging  times,  based  on  the  criterion  selected,  Ergodic  principle, 
etc.;  and  on  the  type  of  data,  analog  vs  sampled  data,  time  signals  or 
power  spectra,  etc.  These  approaches  will  now  he  reviewed. 

In  the  first  approach  (Lumley  and  Panofsky,  1964;  Davenport 
and  Root,  1950;  Bendat  and  Piersol,  1966)  to  questions  of  averaging 
times  the  question  posed  Is  how  long  must  a  time  average  be  in  order 
for  it  to  represent  an  ensemble  average.  The  discussion  considers  a 
random  stationary  analog  signal,  call  it  f(t).  f(t)  could  be  a  direct 
signal  or  some  function  of  a  signal.  The  object  is  to  relate  the  time 
average, 
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and  the  ensemble  average  <f(t)>.  Ensemble  averages  are  denoted  by 
angular  brackets.  The  measure  of  the  difference  between  the  two  is  taken 
to  be  the  enseirble  average  square  difference,  cr2(t)  between  the  two. 


(2)  a2(t) 


f(t+t')  dt' 


<f(t)>]2> 
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a  (t)  is  the  ensenble  variance  of  the  time  average. 

An  expression  for  an  appropriate  averaging  time  is  worked  out 
in  terms  of  B(t'),  the  ensemble  covariance  of  f(t),  where 


(3)  B(t')  =  <[f(t+t')  -  <f(t)>]  [f(t)  -  <f(t)>]> 
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and  the  autocorrelation  p(t)  =  B(t)/B(0).  The  basic  expression  is 
derived  by  first  expressing  Eq.  (2)  first  as  a  double  integral  and 
then  by  putting  it  in  terms  of  the  autocorrelation  to  give 


(4a)  a2 


[f(t+t')  -  <f(t)>]  [f(t+t")  -  <f(t)>]  dt'  dt"> 


(4b) 


■;  t 

|  p(t'-t")  dt*  dt" 
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The  basic  expression  is  simplified  by  transformation  to  sum  and  dif¬ 
ference  coordinates  with  the  subsequent  performance  of  one  of  the  integrals. 
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(5)  t1  =  ,  t2 
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Equation  (6b)  can  be  further  simplified  if  it  Is  the  case  that  the 
autocorrelation  p(ti)  drops  to  zero  in  a  time  much  less  than  T.  In 
that  case  the  second  term  in  the  parenthesis  Inside  the  integral  in 
Eq.  (6b)  Is  negligible.  The  result  expressed  in  terms  of  I,  the 
Integral  scale  of  the  autocorrelation. 
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(8)  T  =  I  , 
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Equation  (8)  is  put  in  a  mo  re  useful  form  by  expressing  a  as  a 
fraction,  e,  of  the  ensemble  mean,  e  =  o/< f>,  giving 


(9),  T  =  2-6101  I 
<f>  e 


,  J  1  .  (  .  , 

Equation  (9)  is  the  basic  result.  It  is  used  to  give  averaging  times 
formulas  for  several  random  quantities  (6en dat  and  Piersol,  1966). 

i  There  are  several  points  worth  indicating  about  this  result.  The 
first  is  that,  it  expresses  an  averaging  time  for  f(t)  in,  terms  of  the 
ensemble  averages  B(0),  <f(t)>  and  ah  integral  of  an  ensemble  average! 
Thtis,  in  order  to  estimate  the  averaging  time  required  for  time  average1 
to  duplicate  the  ensemble  average,  it  is  necessary  to  already  know  < 
some  of  these  ensemble  averages.  To  circumvent  this  problem,1 one  might 
attempt  a  trial  run,  assume  that  the  time  averaged  variance,  mean*  and 
integral  scales  are  approximately  realistic  and  estimate  the  averaging 
time.  If  the  estimated  value  turns  out  to  be  what  has  already  been 
used,  then  the  measurements  are  at  least  consistent  within  the  given 
framework.  i 1 

I  !  |  ,,  1 

A  second  point  is  indeed  more  significant:  Eq.  (9)  says  that  for 
very  long  averaging  times,  the  time  and  ensemble  averages  become 
identical!  ,  Thus,  in  situations  where  it  is  impossible  or  prohibitive 
to  set  up  an  ensemble  of  representative  situations,  then  a  time  average 
can  in  principle  be  made  to  suffice.  It  was  assumed  that  f(t)  was 
stationary,  so  that  there  are  not  any  long  term  trends  to  worry  about! 

1  1  ’  ,  f 

The  final  point  concerns  an  interpretation  of  Eq.  (9).  One  might 
regard  a  long  time  average  as  a  sort  of  ensemble  average  if  the  time 
record  could  be  partitioned  into  successive  segnients  each  uncorrelated. 
Then  each  segment  could  be  regarded  as  an  independent  ensemble  member 
and  the  time  and  ensemble  averages  would  be  comparable.  In  that  sense 
the  integral  scale  of  the  autocorrelation  gives  a  time  proportional  to 
a  decorrelation  time.  The  quantity  i  ,  , 
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might  be  a  little  Wore  representative  of  a  decorrelation  time.  Then 
the  decorrelatibn  time  divided  by  the  square  of  the  desired  precision 
gives  the  requisite  averaging  time. 

Alternatively,  (Bendat  arid  Plersol,  1966)  one  can  state  the  same 
concept  In  different  terms  by  defining  the  equivalent  bandwidth,  W, 
where 

W  =  1/41. 

11  i  ; 

Then  the  tine -bandwidth  product,  WT  is  defined  as  the  number  of  degrees 
of  freedom  and  identified  with  the  number  of  independent  ensemble 
lumbers.  Thus  the  fractional  deviation  from  ensemble  average,  e,  is 
Inversely  proportional  to  the  time- bandwidth  product,  as  might  be  in¬ 
tuitively  expected. 

i  If  one  assumes  a  gaussian  distribution  of  data,  the  equation 


T  = 


_  2  B 


<f> 


can  be  further  reduced.  If  the  quantity  to  be  measured  is  the  mean  square 
value  of  u,  a  gaussian  random  variable,  then  B(0)  =  <(u2-<u2>)2>  * 

2<u2>2,  <f>2  *  <u2>2t  and 


4  1 
~7 


(Lumley  and  Panofsky,  1964).  This  equation,  which  assumes  a  Gaussian 
distribution,  was  applied  to  some  microthermal  data  taken  on  March  23, 
1971  at  the  Rome  Air  Development  Center.  Figure  1  shows  the  auto¬ 
correlation  vs  delay  time  for  a  single  temperature  sensor.  The  auto¬ 
correlation  does  approach  a  zero  value  so  we  assume  that  the  Integral 
scale  does  exist  and  trends  In  data  are  unobservable.  From  this  graph 
we  find  that  I  =  0.1  sec.  and 


T  .  (4)(.l  sec) 

' 


for 


e  =  10% 
e  =  5% 
£  =  2% 


T  =  4C 
T  =  16 


=  40  sec 
,160  secs 
T  =  1000  secs 
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i 
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A  second  approach  (Chamock  and  Robinson,  1957)  to  the  question 
of  averaging  times  considers  sampled  data,  again  asking  how  long  a 
data-taking  duration  is  required.  This  method  considers  a  random 
function  of  time  sampled  N  times  over  a  long  interval  T  from  which  are 
calculated  the  mean 


l  N 

T-Khf< 


and  variance 


V°>  ■  Jr  iI1  <fi  -  f>2 


The  subscript  N  denotes  explicitely  the  number  of  data  samples  used. 

The  N  data  samples  are  divided  into  sections  of  s  samples  each  and  the 
question  asked  is  how  large  must  s  be  in  order  for  the  section  variances, 
averaged  over  the  sections,  to  approximate  the  complete  variance  BN(0). 

The  approach  is  to  consider  the  variances  of  the  individual 
sections.  The  quantity  used  for  comparison  is  B  (0),  the  section 
variances  averaged  over  all  the  sections 


(12)  Bs  (0) 


i  m-1 

1  l 

m  a=0 


,  as+s 

F  l 

s  i=as+l 


,  as+s 

1  l 

5  as+1 


In  Eq.  (12)  there  are  m  sections  of  s=N/m  samples.  The  sub¬ 
script  "a"  indexes  the  individual  sections.  If  the  whole  trace  is  con¬ 
sidered  as  just  one  section  then  s  =  N,  m  =  1  and  a  =  0.  In  that  case 
Bm(0)  is  merely  the  complete  sample  variance  indicated  in  Eq.  (11). 

If  there  is  only  one  sample  in  each  section,  then 


as+s 

s  =  1 ,  m  =  N  and  V  f,  =  f  , 
as+1  J  a  1 

so  that  each  term  in  the  parenthesis  in  Eq.  (12)  is  identically  zero. 
Thus  B-|(0)  =  0. 

The  question  posed  is  what  is  the  shape  of  the  Bc(0)  curve  as 
a  function  of  s?  The  behavior  one  would  like  to  see  is  to  have  the  mean 
of  the  section  variances,  B  (0),  approach  BN(0)  for  s  «  N.  Then  the 


7 


totj1  number  °f  samples,  N,  and  the  averaging  time,  T.,  is  well  above 

fofwhich  The  K3mum  number  of  san*les  would  be  th«t 

ror  which  f(0)  initially,  comes  within  some  preset  region  near  Fw(0) 

F?a  STfUJh10?  !h?  Bsi?^  vs  s  curve  would  be  sinril®»'  to  that  of’ 
rnS»i;,*!f».the  tot?1  number  of  samPles  is  too  small  then  the  Fs(0) 
curve  might  have  a  shape  more  nearly  like  that  of  Fig.  2b.  Thus  it  is 


(a) 


(b) 

Fig.  2.  Variance  versus  averaging  time. 
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the  shape  of  the  1TS{0)  curve  that  determines  the  appropriateness  of 
the  averaging  time.  The  values  of  N  and  s  must  be  sufficiently  large 
with  s  «  N  to  give  the  proper  shape.  In  a  sense  this  is  also  a  self- 
consistency  method  as  described  because  there  is  no  a  priori  way 
stated  for  determining  optimum  values  of  s  or  T. 

Continuing  with  this  approach  (Chamok  and  Robinson,  1957; 
Pasquill,  1962;  Kahn,  1957  and  Smith,  1962)  the  main  points  are 
made  after  Eq.  {12 )_i s  simplified.  To  do  this,  add  and  subtract  the 
long  term  average,  f,  inside  the  parenthesis  in  Eq.  (12),  multiply  out 
the  square  and  simplify. 


,  m-1  ,  as+s  as+s  0 

(13a)  -B  (0)  -1  l  1  l  (frf  -  l  l  (f-f))2 
s  m  a=0  s  as+1  1  s  as+1  1 


,  m-1  ■,  as+s  ,  ,  as+s  ,  as+s  0 

03b)  -MM  C(frf)  -  2(f.-7)  I  l  (f.-f)  +  (I  l  (f.-f))2] 
m  a=0  s  as+1  1  1  s  as+1  1  s  as+1  1 


(13c) 


Equation  (13c)  was  the  starting  point  for  Chamok  and  Robinson. 

In  the  first  term  in  Eq.  (13b)  the  double  summation  is  merely  the 
complete  sum  over  all  samples  so  this  term  is  merely  BN(0).  Thus  Eq.  (12) 
reduces  to  n 


,  m-1  ,  as+s  9 

(14)  *  (0)  *  B  (0)  -  1  l  (I  l  (f.-f))2 
s  N  m  a*0  s  as+1  1 


The  second  term  in  Eq.  (14)  contains  the  difference  between  the  section 
means  and  the  long  term  mean  and  indeed  should  go  to  zero  if  the 
sections  become  sufficiently  long,  leaving  (0 )  as  desired. 

The  second  term  in  Eq.  (14)  can  be  further  simplified  by  writing 
it  in  terms  of  the  covariance  BQ(k),  where 

,  m-1  as+s-k 

<15>  Vk>  ■  srby  J0  ,JS+1  (W7) 
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Thus  multiplying  out  the  terms  in  the  double  summation  in  the  second 
term  in  Eq.  (14)  and  regrouping  using  Eq.  (15)  gives 


(16a)  8S(0)  -  Bn(0)  =■  1  [l-B0(l-s)  +  |  B0(2-s)+...+  §8„(0)  +•••+  B„(s-1)] 


s  so 


S  O' 


s  O’ 


(16b) 


■  r  <si'!  0-44  a  (k). 

s  k=- (s-1 )  s  0 


If  s  and  N  are  sufficiently  large,  then  the  summation  in  Eq.  (16b) 
can  be  approximated  by  an  integral 


(17a)  Bs<0)  -  Bn(0)  -  1  f  (1  -  M>  B(k) 

-s 


dk 


(17b) 


mi  |  o  -Mjf^dk 


assuming  B(k)  is  a  symmetric  function.  The  right  hand  side  of  Eq. 
(17b)  then  has  a  form  identical  with  Eq.  (6b)  if  we  put 


T  k  T„s 

(17c)  -jj-  =  t1  and  -j-  -  T. 


Tq  is  as  before  time  required  for  all  N  samples.  The  result  is 


(17d)  Bt(0) 


bn  (0) 


2  Bn(0) 


t,  B(tJ 

0  ’  r>  p)T  dti 


Equation  (17d)  can  be  derived  in  a  slightly  different  fashion, 
and  a  more  complete  comparison  between  the  approaches  of  Lumley  and 
Panofsky  presented  first  and  of  Chamok  and  Robinson  can  be  presented 
by  interpreting  the  sunr  in  Chamok  and  Robinson's  approach  in  a 
particular  way  and  by  assuming  a  continuous  rather  than  a  sampled 


WWW 


mam#***/*' 


trace.  Thus  consider  the  trace  of  duration  T0  to  be  divided  up  into 
sections  each  of  duration  T..  Then  identify  the  m  intervals  with 
members  of  an  ensenble.  One  then  computes  the  time  average  and  the 
time  variance,  i.e.  the  time  average  square  deviation,  about  the  time 
average  for  each  ensenble  menber.  The  question  asked  then  is  how 
does  this  ensemble  average  of  the  time  variances  depend  on  averaging 
time,  T,  e.g.,  how  long  must  T  be  in  order  for  the  ensemble  average 
time  variance  to  approach  the  ensenble  variance. 

With  the  interpretations  of  the  sums  in  mind  we  can  rewrite 
Eq.  (12)  by  replacing  the  sum  over  “a"  by  an  ensemble  average  and 
replacing  the  sum  over  "i"  by  a  time  average.  Equation  (12)  becomes 


(18a)  Bj(0)  =  <} 


{f(t)  -  t 


f(t)  dtr  dt> 


and  Eq.  (14)  similarly  becomes 


(18b)  Bt(0)  =  B(0)  -  <[} 


f(t)  dt]  > 


The  second  term  on  the  right  hand  side  of  Eq.  (18b)  is  identical  to 
the  exp.assion  in  Eq.  (4a)  and  can  be  similarly  simplified  giving 


(18c)  B^O)  =  8(0)  - 


j  (1  -  4)  P^i)  dti 

o 


which  is  the  same  as  Eq.  (17d),  if  we  identify  Bw(0)  with  the  ensenble 
variance.  This  reinforces  the  concept  of  individual  time  samples 
being  independent  and  being  regarded  as  individual  ensemble  members. 

It  is  interesting  to  note  here  that  the  same  integral  is  used 
in  two  different  situations.  In  Eq.  (6b)  it  gives  the  mean  square 
ensemble  average  difference  between  time  and  ensemble  averages,  denot¬ 
ing  the  approach  of  these  two  averages  with  increased  averaging  time. 
In  Eq.  (18c)  it  gives  the  enseiible  average  of  the  time  variance,  again 
denoting  the  approach  of  the  two  with  increasing  averaging  time. 

There  is  one  other  result  of  Charnok  and  Robinson;  a  procedure 
for  processing  the  function  Bs(0)  in  Eq.  (17b)  and  Fig.  2  to  give 
directly  the  covariance  B(s).  This  expression  is 
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2 

(19a)  B(s)  =  B(0)  -  i  (s2B  (0 ) ) 

3s  5 

2 

This  result  is  obtained  by  differentiating  s  Be(0)  with  b  (0)  taken  from 
Eq.  (17a)  to  give  s  s 


2  3 
(19b)  |_(s2Bs(0))  =  2sBn(0)  -  |  B(k)  dk 

-s 

a2  2 

(19c)  *— j(s  (0))  =  2Bn(0)  -  2B(s) 

dS 


Using  the  fact  that 


(20)  Bn(0)  -  I  |  (frf)2  =  B(0) 

and  rearranging  gives  Eq.  (19a).  Thus,  from  averaging  times  it  is 
possible  to  pull  out  data  on  other  than  just  averaging  times. 

It  is  also  possible  to  express  the  approach  of  the  averaging 
time  to  its  proper  value  as  given  in  Eqs.  (16b)  and  (17a)  using  a 
spectral  description,  (Chamock  and  Robinson,  1957;  Ogura,  1957; 
Pasquill,  1962  and  Smith,  1962).  Thus,  for  example,  introduce  the 
power  spectrum  S(u)  given  by 
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Equation  (23)  indicates  the  well-known  fact  that  a  finite  averaging 
time  acts  like  a  frequency  filter.  The  filter  function  in  this  case  is 
{1  -  [sin(>suiT)/(J5aiT) j2}  which  acts  like  a  high-pass  filter,  passing  fre¬ 
quencies  higher  than  f0  =  1/T.  Thus  the  period  of  recording  should  be 
sufficiently  long  to  faithfully  reproduce  the  lowest  frequencies  desired. 

To  summarize  the  results  so  far,  it  appears  that  an  acceptable 
averaging  time  for  approach  of  time  average  and  variance  to  the 
corresponding  ensemble  averages  might  be  deduced  by  plotting  a  set  of 
data  In  a  manner  indicated  in  Fig.  2a  and  adjusting  the  total  amount  of 
data  and  length  of  sections  until  a  curve  similar  to  Fig.  2a  is  obtained. 
The  time  average  will  then  be  a  good  representation  to  the  ensemble 
average.  One  should  then  check  that  the  section  durations  are  suf¬ 
ficiently  long  that  the  lowest  frequencies  desired  are  reproduced. 

Finally  It  Is  possible  to  obtain  the  covariance  from  the  curve  in 
Fig.  2a. 

A  third  approach  to  averaging  time  considers  measurement  of  power 
spectra  (Blackman  and  Tukey,  1958)  of  random  stationary  analog  data, 
normally  distributed.  Several  steps  of  data  processing  are  used. 

First  the  data  are  processed  to  give  the  covariance  which  is  then 
multiplied  by  a  bellshaped  curve  which  forces  the  product  to  stay  at 
zero  after  some  lag  time  Tm  after  which  the  covariance  has  become  very 
small.  The  product  is  then  transformed  to  give  the  power  spectrum  and 
the  power  spectrum  Is  then  averaged  over  an  ensemble  of  situations.  The 
question  posed  Is  again,  how  long  must  the  original  data  trace  be  in 
order  for  the  power  spectral  value  to  be  within  a  predetermined  range,  in 
dB,  say  ninety  percent  of  the  time. 

The  averaging  time  question  is  answered  quite  simply  for  power 
spectra  because  they  follow  chi-square  statistics.  For  such  a  case 
the  confidence  limits  are  used  to  determine  the  number  of  degrees  of 
freedom,  k,  which  are  In  turn  related  to  the  power  spectral  mean  and 
variance  as  shown  in  Eq.  (25). 


(25)  k  = 


These  spectral  averages  are  in  turn  related  to  the  averaging  time  T  and 
the  cutoff  time  Tm  to  give  the  desired  results.  The  basic  steps  in  the 
derivation  of  this  method  will  be  outlined.  The  reader  is  referred  to 
the  original  source  (Blackman  and  Tukey,  1958)  for  a  more  detailed 
discussion. 

To  continue  with  the  derivation  outline,  then,  start  with  the  time 
trace  X(t)  of  a  single  ensemble  member.  For  that  case  the  covariance 
is  given  by  (neglecting  end  of  trace  effects) 


(26) 


■;/2 

X(t-j)  X(t+J)  dt 

-T/2 


For  this  same 

(27a)  P(<o)  = 

(27b) 


time  trace,  the  power  spectral  density  would  then  be 

oo 

B(t)  D(t)  cos  d 

-00 

-  r/2 

}  j  dx  D(t)  cos  u x  |  X ( t— ^-)  X ( t+j")  dt  . 

-T/2 


D(t)  is  a  bell  shaped  function  with  the  restrictions 

D(t)  =1  x  =  0 

D(t)  =  0  T  -Tm 

where  Tm  is  the  value  of  t  at  which  the  covariance  becomes  neglegibly 
small,  D(t)  eliminate?  noise  arising  from  fluctuations  for  x  >  Tm. 
The  ensemble  average  of  the  power  spectrum  would  then  be 


(28a)  <P(w)> 


<X(t-|)  X(t+|)>dt 


(?Bb) 


dx  D (x )  BQ(x)  cos  wx 


where 


(29)  Bq(x)  =  <X(t+£)  X(t-£)> 


is  the  ensenfcle  covariance  of  X(t). 

The  mean  square  of  the  spectral  density  can  also  be  estimated 
Starting  with  Eq.  (27b),  we  have  for  a  single  trace. 
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00  00 


(30) 


*00  *00 


/2  T/2 

P2(oj)  -  V  f  dT  f  dx'  f  dt  f  dt *  D(x)  D(t')  X(t-£) 
T  J-  J-  -T/2  -T/2 


X(t+J)  X(t'-  J-  )  X(t'+  J-)  cos  WT  COS  UT 1 


Taking  the  ensemble  average  of  both  sides  of  Eq.  (30)  would  give  a 
fourfold  average  under  the  integral  sign.  This  can  be  evaluated 
assuming  the  data,  (the  X's)  are  normally  distributed.  The  result  is 


(31)  <X(t+£)  X(t4)  X(t'+^-)  X(t'-£-)>  =  Bjt-t'+4L4  Bjt-t'-^) 


+  B0(t-f  +  ^-)  B0(t-f  -  ^1)  +  Bq(t)  Bq(t') 


Equations  (30).  and  (31)  are  evaluated  in  the  spectral  domain.  The  result, 
when  combined  with  Eq.  (28a),  is 


(32)  <(P(uj)-<P(w)>)2>  -  Jjrjjy  j  H2(uj,u1)  P2(u1)  du-j 


where  H(w,w|)  Is  a  symmetrized  transform  of  D(t).  Finally  Eqs.  (29) 
and  (32)  are  evaluated  assuming  the  P(u>)  does  not  change  appreciably 
over  a  spectral  range  2ir/Tm.  The  results,  when  inserted  into  Eq.  (25) 
are 


(33) 


2<P(u)= 


<(P(w)-<P(u>)>)  >  2it 


2TC 

>00  £ 
0  P(u)H(w)duD 

2ir 

* 

CP(u))H(a))]2dw 

2T 

7 
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Equation  (33)  Is  the  basic  result.  It  says  that  the  number  of 
degrees  of  freedom  as  determined  from  the  chi-square  distribution  and  the 
desired  precision  is  given  by  2T/Tm.  This  Is  then  interpreted,  as  before, 
as  the  number  of  statistically  Independent  sections  of  data. 

As  another  approach  we  might  indicate  what  one  could  also  apply 
the  chi-square  distribution  to  finding  the  proper  averaging  time  when 
dealing  specifically  with  the  variance  of  a  normally  distributed  stationary 
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random  variable.  For  example,  let  a  time  trace  be  divided  into  sections 

each  Tm  long,  so  that  they  are  essentially  uncorrelated.  Then  the  mean 

variance  formed  from  the  variances  of  the  individual  traces  should 

follow  a  chi-square  distribution.  Thus,  ope  need  merely  use  thq  chi-square 

tables  in  conjunction  with  the  desired  precision  to  determine  the 

requisite  number  of  degrees  of  freedom.  This  then  gives  the  number  of 

sections  of  length  Tm  in  the  total  data  run. 

’  m  i  ■ 

To  sunmarize,  the  stu<fy  of  averaging  times  must  conpem  the 
criterion  for  setting  the  averaging  time,  and  the  exact  quantity 
measured.  Several  examples  appearing  in  the  literature  have  been  con¬ 
sidered  associated  with  the  approach  of  the  time  average  value  of  the 
mean  and  the  variance  to  the  ensemble  value  and  the  precision  of  the 
power  spectrum  and  variance  to  the  desired  values.  1  1 


PHASE  STRUCTURE  FUNCTION  STUDIES 

i  t  I  <  1  ?  ' 

In  the  examination  of  phase  structure  function  various  curves 
have  been  plotted  to  provide  useful  information  for  comparing  results 
and  for  planning  experiments.  There  are  two  types  of  curves.  The 
first  type  shows  spherical  wave  phase  structure  function  as  a  function 
of  separati on  for  various  ranges  ndrmalized  to  refractive  index  parameter, 
Cfa.  This  data  is  intended  for  direct  comparison  with  experimental  results. 
Tne  second  type  of  curve  shows  the  contribution  to  spherical  wave  phase 
structure  funrtion  of  various  portions  of  the  propagation  pathj  These 
curves  of  relative  contribution  versus  distance  are  clotteji  for  various 
separations,  for  various  values  of  the  parameter  2kLg/L  where  x  =  2u/k 
is  the  wavelength,  L  is  the  total  range  and  L0  is  the  turbulence  outer 
scale.  They  are  also  normalized  to  C§.  Also  Included  are  various  paths, 
including  horizontal  and  slightly  inclined  paths  with  the  beam  propagation 
both  upward  and  downward.  These  curves  are  Intended  to  be  useful  in  the 
design  and  interpretation  of  future  experiments  and  in  the  possible  / 
interpretation  of  present  data. '  '  (  1  , 

•  !  1  '  1 
i  i  The  curves  are  based  on  expressions  already  in  the  literature  ,  ' 
(Carlson,  1969)  derived  using  the  ftytov  approximation,  The  basic 
expression  for  D  (p)  is  1  < 


«  i-  i  2 

(34)  Ds(o)  =  8„2  k2  J  [1-0o(kP)]  j  cos 2k~Z^"")2  *N  (r’  dZ  *  d< 


where 


2ir/k  =  light  wavelength 

observation  point  separation  , 

spatial  frequency 

range  1  > 

distance  from  transmitter  along  propagation  path 
refractive  index  spatial  spectrum. 
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I 


l 
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The  Index  spatial  spectrum  used  is  the  Von  Karman  spectrum  given  in 


The  index 
Eq.  (35). 


(35a)  *„(«) 


0.033  Cjj(H) 

(k2  ♦  1/Lq) 


,  -(Klt0/5.92)i 

OH  e  0 


(35b) 


-  0.033  cj^(H)  (<2  +  l;2)-’V6 


where  i  -  turbulence  inner  scale,  H  =  height. 

i  ' 

! 

I*  .  - 

The  approximate  form  in  Eq.  (35b)  was  used  because  the  phase  structure 
function  Is  not  generally  measured  at  the  separations  comparable  with 
the  inner  scale,  and  the  curves  are  for  comparison  with  measured  values. 

For  the  slaht  paths  the  structure  parameter  was  taken  to  vary  with  height 
to  the  power  (-4/3),  (Wyngard  et  al,  1971)  and  the  outer  scale  proportional 
to  altitude.  These  are  valid  In  the  first  500  ft  for  unstable  air  and  well 
developed  turbulence. 

(36a)  C'j(H)  =  cg(H0)(H/Ho)'4/3 


(36b)  Lq(H)  «  L0(H0)({j-) 

H  7  H,-H0 

<*'>  .fc-1 

O  O 


i  where 


Hq  *  transmitter  height 
Hl  =  receiver,  height. 


Equation  (36c)  written  with  distance  along  the  path  as  a  parameter  was 
the  actual  expression  used.1 

The  calculations  were  performed  on  a  digital  computer  using 
expressions  derived  from  Eqs.  (34)  -  (36).  To  find  the  Incremental 
pathlength  contributions  to  Ds(p)  the  order  of  integrations  in  Eq.  (34) 
was  reversed,  the  resultant  Inner  integral  giving  the  desired  Incremental 
path  contributions.  The  calculations  were  made  using  normalized 
dimensionless  parameters 


i 


i 


(37a)  u  =  L0(H0)k 

(37b)  v  =  Z/L 

(37c)  y  2k  (HQ)/L 

(37d)  CQ  =  C2(HQ)  k7/6L11/6  =  C*(0)/0.13 


(Tatarski,  1961)  giving  to  Eq.  (34)  the  form 


1 

(38a)  Ds(p/LQ)  =  {0.033  x  8*2C0  x  (^-)t>/6}  x  j  dv  F(v) 

o 


where 

(38b) 


F(v) 


du(1-J  (P^)) 
o 


X  (1  + 


H,  -  H. 


(- 


H. 


'-)) 


-4/3 


Of  the  reduced  parameters,  u  is  a  reduced  spatial  frequency,  v  is 
reduced  range,  y  is  the  Fresnel  number  for  an  aperture  the  size  of  the 
outer  scale,  and  C0  is  proportional  to  the  spherical  wave  log  amplitude 
variance,  often  used  as  a  scaling  parameter. 


Figures  3a  and  3b  show  the  contribution  to  the  normalized  phase 
structure  function  for  8  separations  (normalized)  and  for  y  =  10+3.  The 
normalized  function,  F(v)  is  insensitive  to  the  parameter  y  and  hence 
Fig.  3a  covers  all  of  the  experimentally  interesting  range  of  y.  Fig. 

3b  shows  that  y  has  a  larger  (but  still  small)  effect  on  F(v)  for  the 
smaller  separations  (p/L0).  This  means  that  the  range  dependence  in 
F(v)  Is  very  nearly  independent  of  wavelength,  the  values  1()3  ^y  =  <  io5 
cover  both  0.6p  and  10p  for  rangas  of  Interest.  The  wavelength  dependence 
is  in  the  y5/6  part  of  the  normalization  constant  in  front  of  the  Integrals. 
From  Figs.  3a  and  3b  It  is  evident  that  the  most  Important  contributions 
to  Ds  are  made  near  the  receiver  for  a  horizontal  path. 
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Differential  spherical  wave  phase  structure  function  contribution 
vs  normalized  range  for  five  normalized  separations. 


Fig.  3b.  Differential  spherical  wave  phase  structure  function  contributi 
vs  normalized  range  for  three  normalized  separations. 


Figures  4  and  5  show  the  normalized  Ds  integrand  for  an  upward  and 
downward  slant  path  respectively.  In  the  upward  path  the  curve  amplitude 
has  "saturated",  so  that  a  more  nearly  constant  contribution  to  range 
results.  In  the  downward  path  the  turbulence  effects  are  strongly 
concentrated  at  the  receiver. 

Having  the  functions  in  Figs.  3a  and  3b  it  is  a  simple  matter  to 
calculate  the  complete  phase  structure  function. 

The  unnormalized  phase  structure  function,  for  direct  comparison 
with  experimental  results,  is  presented  in  Figs.  6a,  6b,  and  6c  for  the 
three  ranges  of  immediate  interest. 

Having  developed  the  programs  to  calculate  range  dependence  of 
the  phase  structure  function,  it  was  also  possible  with  only  minor 
modifications  to  calculate  the  range  dependence  of  the  small  aperture 
angle  of  arrival  correlation  functions. 

The  small  aperture  elevation  (IL)  and  the  azimuth  (Bg)  angle  of 
arrival  correlation  functions  are  related  to  the  phase  structure  function 
by 


(39) 


B  (p)  =  \ 
°  2kZ 


320s(p) 


(40) 


V»> 


aOs(p) 

3p 


Equations  34  ,  39,  and  40  may  be  used  to  find  the  contribution 
to  Ba  and  Bfi  from  various  ranges,  and  (for  simplicity)  a  horizontal 
path. 


(41) 


Bg(^)  =  4ir2 


1 

x  .033  C2(Ho)  -^/3  |  G(v)  dv 
Lo  o 


(44-i})(1+4)’11/6  du 

3  v 
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Differential  spherical  wave  phase  structure  function  contribution 
vs  normalized  range  for  three  normalized  separations. 


o 

o 


:5 


Fig.  6b.  (Cont. ) 


KV4KM 


•xrnmum^p 


(42)  Bo(j^)  =  4tt2  x  .033  C2(Hq) 


d 


1}) 


x  (1 


+ 4ri,/6 


du  -  G(v). 


Typical  examples  are  shown  In  Figs.  7a  and  7b. 

(17  s i gni fl HaSt 'figj res I^Hthmetir^A^6  "w  r°uble  prects1on 
integration  subroutine  was  used^to  do  if  P°lnt-G*ussian  quadrature 
numerical  Integration  11  Sits  iJe^uS  the  Inner  integrations.  The 
accuracy.  minis  are  such  as  to  yield  better  than  1  % 
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BIBLIOGRAPHY  OF  OPTICAL  PROPAGATION 
IN  A  TURBULENT  ATMOSPHERE 


The  final  item  is  a  continuation  of  the  bibliography  included 
in  Ohio  State  University  Report  RF  2880-2,  and  includes  journal  articles 
dealing  with  the  effects  of  a  turbulent  atmosphere  on  electromagnetic 
wave  (generally  light  beam)  propagation.  These  articles  are  listed  by 
year  and  then  alphabetically  by  author.  One  of  the  major  features  of 
this  revised  bibliography  is  the  inclusion  of  Russian  articles  in  Izv. 
Vuz.  Radio  Fizika  that  have  not  been  translated  into  the  English  edition 
Soviet  Radio  Physics.  These  original  Russian  articles  have  English  titl 
and  abstracts. 

There  was  also  some  additional  information  on  authors  and  journals 
compiled  from  the  complete  bibliography.  This  information  is  presented 
in  two  additional  listings. 

List  1  tabulates  the  journals  surveyed,  the  years  and  the  number 
of  entries  found  for  the  complete  bibliography. 

List  2  gives  an  alphabetical  tabulation  of  the  first  authors  in 
the  complete  bibliography  and  the  years  that  they  published. 
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SUMMARY 


In  this  interim  report  three  areas  of  interst  have  been  considered: 
averaging  times  for  random  data,  phase  structure  function  data  and  a 
continuation  of  a  bibliography.  The  averaging  time  study  reviewed 
approaches  contained  in  the  literature  showing  that  the  averaging  time 
depends  on  the  quantity  of  interest  and  on  the  particular  criterion.  The 
phase  structure  function  computations  gave  curves  intended  to  be  useful 
in  experiment  design  ard  data  interpretation.  The  bibliography  presented 
a  listing  of  recent  journal  articles  and  author  and  journal  lists. 
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